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ABSTRACT 


In  this  report,  the  ability  of  o  radar  to  resolve  overlapping  signals  is  exotnined.  It 
is  shovm  that  the  primary  limitation  on  resolvobility  is  the  unpredictable  difference 
between  the  actual  received  signal  and  the  waveform  for  which  the  target  receiver 
is  matched  (the  distortion  in  the  received  signal).  The  effects  of  signal  distortion  on 
the  shape  of  the  ‘ambiguity  function*  are  examined,  based  upon  a  wideband  random 
process  model  for  the  distortion.  The  statistics  of  the  ambiguity  function  side-lobe 
region,  in  the  presence  of  the  distortion,  are  related  in  a  simple  manner  to  the  sta¬ 
tistics  of  the  distortion  process.  The  optimum  two-target  resolver  for  the  distorted 
signal  is  then  derived  and  its  performance  is  examined  for  two  different  cortditions 
of  o  priori  knowledge  concerning  the  signal:  (1)  All  parometers  of  both  target  returns 
krsown  exactly,  with  only  the  presence  or  absence  of  the  second  target  uncertoin, 
and  (2)  the  phase  and  amplitude  of  either  or  both  returns  unknown,  butall  other  pa¬ 
rameters  exactly  known.  These  two  cases  moke  evident  the  comparative  seriousness 
of  wide-band  distortion  and  lock  of  o  priori  knowledge.  Probability  of  detection  vs 
probability  of  folse-alarm  curves  are  derived  for  both  eases. 

The  performance  of  two  nonoptimum  two-torget  resolvers  which  ore  simpler  from  a 
circuitry  viewpoint  than  the  optimum  resolvers  is  examined.  Included  in  the  non- 
optlmum  resolver  is  the  usual  ‘matched*  (to  a  single,  isolated  return)  filter.  It  is 
quantitatively  demonstrated  that  for  nearly  complete  overlap  In  time  of  the  received 
waveforms  from  the  two  targets  (even  though  for  ‘compressed*  signals  the  targets 
may  be  separated  by  more  than  the  width  of  the  main  lobe  of  the  ambiguity  function) 
the  resolvability  of  a  small  target  in  the  presence  of  a  large  target  tends  to  become 
independent  of  the  radar  sensitivity.  Thus,  for  sufficiently  large  cross-section  dif¬ 
ference,  resolution  con  be  achieved  only  by  lowering  the  ambiguity  function  side- 
lobe  uncertainties,  os  one  would  intuitively  expect.  The  optimum  vs  nonoptimum 
resolver  performonces  differ  primarily  in  the  rote  at  which  targets  become  resolvable 
with  decreasing  overlap. 
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A  QUANTITATIVE  EXAMINATION  OF  THE  RADAR  RESOLUTION  PROBLEM 


I.  INTRODUCTION* 

A.  The  Radar  Resolution  Problem 

Perhaps  the  most  vexing  problem  for  the  radar  designer  is  obtaining  adequate  resolution  (the 
detection  of  targets  in  the  presence  of  other  reflecting  objects).  For  operation  at  relatively  short 
ranges,  the  principal  resolution  problem  is  usually  resolving  targets  in  a  background  of  more  or 
less  continuous  ground  or  sea  clutter.  At  longer  ranges,  the  earth's  curvature  obscures  most 
of  the  sources  of  clutter  return,  and  then  the  resolution  problem  is  to  resolve  discrete  tar¬ 
gets  from  one  another,  particularly  at  frequencies  above  UHF  at  which  aurora  and  meteors  are 
not  important.  This  report  represents  an  investigation  of  the  discrete  target  problem. 

Impetus  for  this  investigation  arises  from  consideration  of  several  contemporary  radar 
problems  which  demand  the  resolution  of  individual  targets  of  small  radar  cross  sections  at 
ranges  of  hundreds  of  miles.  Such  problems  seem  inevitably  to  be  complicated  further  by  the 
fact  that  the  target  of  interest  is  much  smaller  in  radar  cross  section  than  the  interfering  target.t 
Even  a  rudimentary  consideration  of  this  kind  of  resolution  problem  makes  it  clear  that  the 
spread  of  anticipated  cross  sections  determines  the  feasibility  of  different  approaches  to  the  prob¬ 
lem.  If  very  great  differences  (several  tens  of  db's)  in  cross  section  must  be  resolved,  the  de¬ 
signer  has  little  choice.  The  only  radar  "dimension"  which  can  be  made  entirely  free  of  "side 
responses"  is  that  of  range,  since  both  the  Doppler  response  and  the  angular  response  of  a  radar 
possess  "side  lobes,"  and  the  side  lobes  of  the  radar  response  to  a  large  target  may  obscure  the 
main  response  to  a  small  target.  The  straightforward  method  of  achieving  this  "absolute"  reso¬ 
lution  in  range  is  to  use  a  signal  consisting  of  a  single  short  pulse.  However,  the  ability  of  a  ra¬ 
dar  to  detect  the  presence  of  a  target  and  to  estimate  target  parameters  in  a  background  of  white 
thermal  noise  is  determined  by  the  ratio  of  the  energy  E  received  from  the  target  during  one 
decision-making  interval  to  the  noise  power  per  cycle  N^.  Obviously,  achieving  sufficient  energy 
in  very  short  pulses  can  result  in  severe  demands  on  radar  components. 

Furthermore,  it  is  by  no  means  clear  that  for  all  radar  applications,  range  is  the  most  useful 
dimension  in  which  to  achieve  good  resolution.  For  example,  consider  the  problem  of  resolving 
newly  injected  satellites  from  their  boosters.  The  impulsive  separating  mechanism  causes  the 


*Thi(  report  it  a  rovition  of  lh«  thotit  tufamittod  by  iho  author  to  tho  Doportmont  of  Electrical  Engineering,  M.I.T., 
on  19  May  1962,  in  partial  fulfillment  of  the  requirements  for  the  degree  of  Matter  of  Science. 

t  An  example  of  a  problem  of  this  cion  it  that  of  obtaining  early  resolution  of  o  newly  lourtched  satellite  from  its 
booster  rocket  to  enable  precise  prediction  of  the  satellite  oibit  from  the  radar  measurement. 
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two  bodies  to  differ  instantaneously  in  velocity  upon  separation,  while  they  gradually  separate  in 
range.  In  principle,  then,  the  bodies  are  potentially  resolvable  in  velocity  (Doppler)  before  they 
are  resolvable  in  range. 

Thus,  it  is  often  impractical  and/or  undesirable  to  attempt  to  achieve  high  resolution  only 
in  range  by  use  of  short  pulses  exclusively.  The  question  which  naturally  arises  is  how  much 
"limited"  (by  side-lobe  responses)  resolution  can  be  obtained  in  the  dimensions  in  which  such 
side  lobes  must  exist,  such  as  Doppler  frequency  or  angular  dimensions.  Once  the  possibility 
of  utilizing  "limited  resolution"  is  raised,  it  is  also  appropriate  to  examine  the  utility  of  so-called 
"pulse-compression"  techniques,’’'  which  create  range  side  lobes  in  addition  to  those  in  the  other 
dimensions. 

Therefore,  we  shall  investigate  here  fundamental  and  practical  limitations  upon  the  degree 
of  "limited"  resolution  obtainable  from  a  pulsed  radar,  with  emphasis  on  those  using  pulse- 
compression  waveforms. 

B.  The  Direction  and  Scope  of  this  Study 

We  shall  describe  two  aspects  of  the  limited  resolution  problem.  First,  we  will  try  to  de¬ 
termine  the  fundamental  limitation  on  resolution  of  returns  from  targets  that  are  sufficiently 
"close"  (in  terms  of  the  observable  radar  coordinates,  such  as  range,  Doppler  and  angular  meas¬ 
urement)  so  that  their  return  signals  "overlap."  In  particular,  we  will  derive,  through  a  simple 
but  informative  model,  the  probabilities  that  we  will  correctly  or  incorrectly  announce  the  pres¬ 
ence  of  a  second  target  in  the  side- lobe  structure  of  a  first. 

Investigation  into  fundamental  resolvability  specifies  an  optimum  circuit  for  resolving  the 
targets.  However,  even  in  the  happily  abstract  environment  of  this  study,  the  hardware  realiza¬ 
tion  of  the  optimum  resolver  must  be  admitted  to  be  highly  complex.  Thus,  we  will  also  examine 
the  degree  to  which  certain  less  complicated  processing  techniques  compare  to  the  optimum  re¬ 
solver  in  performance.  Included  among  the  processors  studied  is  the  usual  "matched"  filter,  a 
filter  with  impulse  response  which  is  the  undistorted  return  from  a  single,  isolated  target  re¬ 
versed  in  time. 

C.  The  "Ambiguity  Function”  and  its  Relationship  to  Resolution 

Prior  to  the  early  1950*8,  investigations  of  the  resolving  capabilities  of  radar  appear  to  have 
been  entirely  qualitative,  and  almost  exclusively  directed  at  the  generation  and  reception  of 
shorter  pulses.  For  example,  examination  of  the  index  to  the  M.I.T.  Radiation  Laboratory  Series 
reveals  only  a  single  reference  to  target  resolution:  a  brief  statement  about  the  importance  of 
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range  resolution  in  airborne  radar  for  high-speed  aircraft. 

Woodward  appears  to  have  first  explored  the  more  general  problem  of  resolution  in  more 
than  one  dimension,  discussing  resolution  in  both  range  and  Dc^pler.  He  is  generally  credited^ 
with  first  pointing  out  that  a  fundamental  measure  of  the  ability  of  a  radar  to  resolve  signals 
from  two  targets  with  different  parameters  (specifically  range  and  Doppler,  in  Woodward's 
treatment)  is  the  shape  of  the  cross-correlation  function  of  the  returns  from  two  targets  with 
differing  parameters. 

*A  brief  deieription  of  Mmo  typical  wavaforms  for  these  techniques  is  given  in  Sec.l-A;  for  mere  information, 
the  reader  is  referred  to  the  BiUiogrophy  (Refs.  1, 2,  3).  Such  waveforms  are  characterized  by  a  time  dimthm- 
bondwidth  product  much  larger  than  unity  (TW  »  1,  where  W  is  In  cycles  per  second,  T  in  seconds). 

tSec.7.1  ofRef.5. 
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To  formally  define  this  function,  we  assume  that  the  return  from  a  target  is  representable 
by  the  real  part  of  some  complex  waveform,^  p(t,  where  generality  can  be  maintained  con¬ 
veniently  by  regarding  7^  as  an  abstract  "vector,"  the  components  of  which  represent  the  effect 
of  the  target  parameters  of  Target  Number  1  on  the  signal.  For  this  notation  the  cross¬ 
correlation  function^  between  two  targets  is  defined  by 


(1) 


where  the  asterisk  denotes  the  complex  conjugate,  and  the  p's  are  generally  normalized  so  that 

xcyi-yj)  =  1: 

j*  |p(t,7i)l^  dt  =  1  .  (2) 


This  function  is  known  as  "Woodward's  Ambiguity  Function,"  or  "Woodward's  Uncertainty 
Function." 

If  we  let 

=  p(t)  .  (3) 


for  targets  differing  only  in  range  and  Doppler  and  for  narrow-band  waveforms, 
imate^  pit,  ^2) 


pit,^^)  =  -  t)  e 


we  can  approx- 


(4) 


where  t  is  related  to  the  range  difference  R  by  t  =  2R/c,  and  Wp  represents  the  usual  narrow- 
band  approximation  to  the  Doppler  shift.  Then  the  range-Doppler  form  of  the  ambiguity  function 
is 

X(t,  Wp)=J  p{t)p*(t-T)e  dt  .  (5) 


As  Woodward  indicates,  this  function  plays  an  important  role  in  resolution  questions,  because 
its  amplitude  expresses  the  degree  to  which  signals  from  different  targets  differ  in  structure  in 
the  dimensions  included  in  the  7  vectors.  When  |x|  is  near  unity,  the  received  signals  are 
nearly  identical  in  the  dimensions  concerned,  and  one  would  anticipate  that  when  the  two  signals 
are  immersed  in  a  noise  background,  distinction  would  be  difficult.  Conversely,  a  value  of  |  X  | 
near  zero  indicates  distinctive  signals,  and  indicates  less  difficulty  in  resolution.^ 


tComplex  functions,  dsnotsd  primarily  by  Greek  symbols,  will  be  used  throughout  this  report,  since  the  mathemat¬ 
ical  natation  is  considerably  simplified.  Vvhen  we  speak  ^  a  function  a(t)  we  imply  a  function  z(t)  «  x(t)  +  iy(t), 
whose  real  port  x(t)  Is  the  actual  process  of  concern,  and  whose  imaginary  part  y(t)  is  the  Hilbert  tronsform  o  fx(t). 
This  particular  specification  for  y(t)  renden  the  spectrum  of  z(t)  single  sided  (Z(w)  0,  u  <  0],  and  allows  one  to 

write  <(t)  modulated  by  a  CW  carrier  of  angular  frequency  uo  simply  os  z(t)  exp  (juof )•  Other  salient  features  of 
this  notation  will  be  discussed  os  required;  the  unfamiliar  reader  is  referred  to  Woodward,'  Sec.  2. 9;  Kelly  and 
Reed,6Doob.7  ^ 

t  The  use  of  the  symbol  X  is  by  no  means  universal;  many  symbolisms  ore  in  use.  This 

notation  follows  Helstrom.  vr°r 

§  The  Doppler  shift  is  actually  a  spectrum  "stretch,"  rather  than  a  spectrum  shift.  The  latter,  however,  is  gen¬ 
erally  corrceded  a  sotlsfectory  approximotion  for  narrow-bond  signals. 

HAIthough,  os  we  shall  see  for  |  X|, the  Inevitable  urKertainty  in  |  X|  due  to  sigtwl  distortion  ond  system  errors 
becomes  more  of  a  fundamental  limitation  than  |  X|  itself. 
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Further  discussion  of  the  role  of  X  in  range  and  Dc^pler  resolution  is  given  in  Woodward^ 

and  Siebert  Kelly*^  has  explicitly  extended  the  ambiguity  function  to  the  dimensions  of 
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range,  radial  velocity  and  radial  acceleration,  and  Mayo,  st  &l->  have  investigated  its  prc^er- 
ties  as  a  function  of  range,  Doppler  and  two  orthogonal  angular  measurements,  thus  including 
certain  antenna  effects  which  become  important  with  extremely  wide-band  signals. 

D.  Further  Studies  of  the  Resolution  Problem 

In  the  previously  mentioned  studies  of  radar  resolution  no  attempts  were  made  to  quantita¬ 
tively  assess  the  resolution  capability  of  radar  (probabilities  of  detecting  a  single  target  in  the 

g  q 

presence  of  others).  Helstrom  '  was  the  first  to  attack  the  problem  quantitatively,  the  later 
work  was  more  nearly  a  model  of  the  usual  radar  situation.  He  considered  the  form  and  per¬ 
formance  of  the  maximum -likelihood  detector  for  the  reception  of  two  signals  of  known  separation 
in  time  and  frequency,  but  with  individual  phases  and  amplitudes  unknown,  and  the  time  of  ar¬ 
rival  of  the  pair  unknown.  The  two  signals  can  be  written  (in  our  notation)  with  Signal  1  expressed 
as 

«l(t)  =^l^(t-y 

and  Signal  2  as 

^^(t)  =^2P(t-t^-T)  e  °  ° 

where  and  complex  amplitudes,  and  and  t^  are  assumed  unknown.  Helstrom 

concludes  that  as  |  ->  «,  the  detection  Signal  2  in  the  presence  of  Signal  1  becomes  insen¬ 
sitive  to  1^^  I,  and  the  cost  in  required  increase  in  returned  energy  (from  Target  2)  for  a  given 
Pjj  and  Pp  is  only  1  or  2db  for  |x|^  <  0.1,  even  though  the  amplitude  of  the  interference  of  Sig¬ 
nal  1  at  the  location  of  Signal  2  jX|  is  essentially  infinite.  This  result  seems  rather  contrary 
to  experience.  As  will  be  proven  later,  it  is  a  consequence  of  assuming  the  form  of  the  signals 
|i(t)  and  their  separation  to  be  exactly  known  -  a  situation  which  does  not  prevail  in  practice. 

Fowle,  Kelly  and  Sheehan^  ^  discuss  a  situation  in  which  there  are  many  targets  whose  am  • 
biguity  functions  overlap  in  a  manyfold  manner.  By  assuming  that  the  multitarget  background 
resembles  Gaussian  noise,  an  expression  for  the  maximum  number  of  targets  resolvable  (on  the 
average)  was  derived,  in  terms  of  required  signal-to-noise  ratio  and  the  cross-section  spread 
present. 

In  a  recent  paper,  Nilsson^^  qualitatively  considered  the  optimum  resolving  system  for  a 

Bayes  mean-square  error  criterion  for  a  situation  in  which  there  is  a  known  number  n  of  point 

targets  which  differ  only  in  range.  He  also  extended  his  results  to  include  an  unknown  number 

of  targets,  adopting  a  general  penalty  function  a^^  for  guessing  i  targets  when  j  are  present.  He 

makes  no  calculations  of  the  optimum  processor's  performance,  nor  does  he  attempt  to  extend 
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results  to  resolution  in  dimensions  other  than  range.  Preston  has  also  briefly  considered  a 
"near-optimum*  resolver,  but  again  in  only  qualitative  terms.  Thus,  we  can  see  that  the  pub¬ 
lished  quantitative  work  on  this  radar  resolution  problem  is  quite  meager.'*' 


*  Sines  the  original  drafting  of  thb  rsport,  W.L.  Root  ho*  alto  publiihod  a  quantitative  examination  of  the  reso¬ 
lution  problem  in  "Radar  Resolution  of  Closely  Spaced  Targets,"  Trans.  IRE  PGMIL,  MIL-d  1^  (April  1962). 
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n.  THE  PROBLEM  MODEL 


In  this  chapter,  the  considerations  that  led  to  the  model  of  the  resolution  problem  that  was 
studied  are  reviewed.  First,  the  many  types  of  waveforms  In  use  or  contemplated  are  briefly 
reviewed  and  their  common  features  Indicated.  The  many  possible  sources  of  distortion  that 
exist  for  such  radar  signals  are  pointed  out,  and  a  statistical  model  of  their  aggregate  effect 
proposed.  This  signal  distortion  is  then  related  to  the  resulting  ambiguity  fimctlon  uncertainty, 
and  some  simplifying  assumptions  about  the  problem  model  found  feasible. 

A.  Some  Characteristics  of  Common  Radar  Waveforms 

The  details  of  modulation  techniques  used  or  prc^osed  for  radar  differ  considerably,  but  the 

more  widely  used  techniques  can  be  categorized  approximately. 

5  11 

Since  we  can  establish  ’  that  the  principal  lobe  of  the  range-Doppler  ambiguity  function  is 
roughly  of  an  extent  in  time  equal  to  the  reciprocal  of  the  bandwidth’*'  l/W  and  of  an  extent  in 
Dqspler  frequency  roughly  equal  to  the  reciprocal  of  the  signal  duration  l/T,  we  see  that  the 
so-called  "duration-bandwidth"  product  TW  is  significant.  For  example,  TW  approximately 
specifies  the  ratio  of  the  transmission  duration  to  the  extent  in  time  of  the  main  lobe  of  the  am¬ 
biguity  function  (the  "pulse-compression  ratio"). 

We  will  concentrate  our  attention  on  signals  of  large  TW  product.  It  is  convenient  to  classify 
such  signals  by  whether  the  bandwidth  is  achieved  by  modulation  within  one  long  continuous  en¬ 
velope,  or  by  using  a  number  of  individual  pulses.  The  former  class  will  be  referred  to  as 
"coded-pulse"  signals;  the  latter  as  "pulse-train"  signals.  Examples  of  coded-pulse  signals  are: 
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(1)  Pulses  with  frequency  modulation  '  within  the  pulse,  as  indicated 
in  Fig.  1(a),  with  achievable  TW  products  limit^  only  by  hardware 
constraints. 

(2)  The  so-called  "phase  codes,"  which  make  use  of  discretely  switched 
phase  within  the  pulse  envelope,  of  which  perhaps  the  best  known  is 
the  "Barker*!^  code  of  order  13,  as  indicated  in  Fig.  1(b),  with 

TW  a  13.  Greater  TW  products  are  possible,  of  which  t)ie  so-called 
"shift-register"  codes  are  an  example. 

Typical  of  pulse-train  signals  are  those  which  are: 

(1)  Regularly  spaced,  equal-width  pulses,  as  in  Fig.  l(c),^^  or 

19 

(2)  Trains  with  nonuniform  spacing,  width,  etc. 

We  might  vary  any  of  the  potential  train  parameters  in  "pulse  train"  signals,  but  hardware  con¬ 
siderations  usually  lead  to  varying  only  the  interpulse  spacing  as  in  Fig.  1(d). 

Combinations  of  the  various  types  may  also  have  desirable  prc^erties  for  some  applications 
as  well,  of  course,  but  will  not  be  explicit^  considered  here. 

Note  that  the  signals  of  these  classes  used  in  radar  are  of  a  binary  nature;  the  amplitude  is 
either  zero  or  a  constant  nonzero  value,  due  to  the  difficulty  in  making  efficient  radar  transmit¬ 
ters  which  operate  as  linear  amplifiers  (high  efficiency  requires  operation  in  a  saturated  and 
therefore  constant  amplitude  mode).  This  fact  is  of  analytic  utility  and  will  be  used  in  the  anal¬ 
ysis  to  follow. 


*W  it  in  eyelet  per  teeond  throughout  thit  report. 
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B.  DlBtortlon  of  Radar  Signals  -  Causes  and  Effects 

Helstrom's^  and  Nilsson results  indicate  that  if  one  knows  the  form  of  the  return  signal 
(and,  consequently,  the  shape  of  the  ambiguity  function)  exactly,  the  effects  of  the  ambiguity 
function  side  lobes  of  a  large  target  can  be  completely  removed  by  noting  the  target  "position" 
and  subtracting  an  appropriately  positioned  and  weighted  ambiguity  function.  However,  any  un¬ 
certainty  in  the  form  of  the  received  signal,  resulting  in  an  uncertainty  in  the  shape  of  the  am¬ 
biguity  function,’^  partly  nullifies  this  cancellation.  Furthermore,  since  the  size  of  the  effect  of 
the  target  is  given  by  the  product  of  its  cross-section  magnitude  and  ambiguity  function  uncer¬ 
tainty  (recall  that  the  ambiguity  function  is  normalized),  no  matter  how  small  the  uncertainty  in 
the  ambiguity  function,  a  large  enough  cross  section  exists  to  render  its  effect  significant. 

Thus,  we  must  be  concerned  with  the  effects  of  even  very  small  ambiguity  function  uncertain¬ 
ties  (values  of  perhaps  -30db  or  so  being  "state  of  the  art"),  and  consequently,  the  possibility 
of  even  small  distortions  of  the  signal  returned  to  the  radar  as  well. 


Fig.  2.  Conceptual  diagram  of  a  typicol  radar  (ituotion. 

The  possible  sources  of  signal  distortion  are  many.  To  appreciate  this  fact,  we  need  only 
consider  the  route  of  a  radar  signal  from  the  generation  of  the  modulation,  to  the  target,  and 
back  to  the  radar  receiver,  as  indicated  in  Pig.  Z.  With  no  attempt  at  being  exhaustive,  we  can 
enumerate  some  typical  distortion  sources. 

(1)  Both  the  phase  and  amplitude  of  the  transfer  coefficients  of  the  trans¬ 
mitter  and  receiver  may  not  be  linear  or  flat,  respectively.  The  form 
of  the  distortion  will,  furthermore,  have  a  time-varying  component 
due  to  such  causes  as  thermal  effects  and  power  supply  variations. 

(2)  The  frequency  characteristics  of  the  transmitting  and  receiving  anten¬ 
nas  as  "seen"  by  the  target  depend  not  only  upon  the  feed  network 
characteristics,  but  also  PPbn  the  angle  of  the  target  from  the  pointing 
direction  of  the  antenna,^”’*'  due  to  the  finite  propagation  time  of  en¬ 
ergy  across  the  antenna  from  angles  off  broadside.  Further,  the 
characteristics  will  depend  upon  the  antenna  pointing  angle,  relative  to 
some  fixed  reference,  due  to  such  causes  as  rotary  joint  characteristics, 


*  Note  that  such  "uncertainty"  may  be  either  truly  a  matter  of  ignorance,  or  merely  a  matter  of  convenient  neglect 
of  available  information. 


7 


reflections  from  nearby  objects,  and  with  phased  array  antennas, 
the  so-called  "build-up"  phenomena  (Part  3,  Ch.  2  of  Ref.  21),  due 
to  the  finite  time  of  wavefront  propagation  across  the  array. 

(3)  The  effects  of  the  target  on  the  radar  waveform  are  many  and  varied. 

First,  the  Doppler  shift  introduced  by  the  component  of  the  target 
motion  relative  to  the  radar  will  actually  be  a  spectrum  "stretch,* 
rather  than  the  spectrum  displacement  usually  assumed;  the  re¬ 
flected  frequency  oij.  is  related  to  the  instantaneous  incident  frequency 
and  the  ratio  of  the  velocity  component  in  the  radar  direction* 

Vj,  to  the  velocity  of  propagation  c  1^ 

1  -v^c 
"r  =  “i  TT^ 

or,  to  a  first  order, 

r  1 

“r  =  “il^-— I  •  (6) 

The  simple  constant  displacement  frequently  used  for  analytic  purposes 
is  only  accurate  for  very  narrow -band  signals. 

«r  =  “i  “d  * 


where  Wq  is  assumed  to  be  a  constant  related  only  to  the  carrier 
frequency 


Furthermore,  Vr/c  can  be  taken  as  constant  only  if  the  transmission 
time  is  very  short  compared  to  any  target  acceleration.  Finally,  if  the 
time  extent  of  the  target  is  of  a  dimension  comparable  to  the  reciprocal 
of  the  signal  bandwidth  and  nonsymmetrical  as  viewed  from  the  radar, 
the  target  cross  section  may  vary  in  both  phase  and  amplitude  with  fre¬ 
quency  over  the  signal  bandwidth. 

Not  only  are  there  many  sources  of  distortion  in  the  path  of  a  radar  signal,  but  some  of  these 
sources  are  capable  of  producing  distortion  of  a  highly  structured,  wide-band  nature.  If  the  sig¬ 
nal  is  narrow-band,  the  receiving  filters  will  remove  most  of  the  distortion  structure  from  the 
signal.  However,  for  resolution  at  long  ranges,  we  are  interested  in  wide-band  signals,  and  must 
consequently  be  prepared  for  similar  distortion. 

An  attempt  to  explicitly  describe  the  effects  of  the  many  possible  distortions  is  difficult. 
Furthermore,  since  many  of  the  sources  of  signal  distortion  are  of  a  random  nature,  the  use  of 
a  statistical  model  suggests  itself.  The  fact  that  we  expect  a  time  structure  to  the  distortion, 
leads  us  to  consider  the  distorted  signal  from  the  i^^  target  to  be  of  the  form 


^Kt,y.)  =  (3j  +  «j(t))p(t,yj)  ,  (8) 

where  8^  is  the  complex  target  cross  section  and  8.(t)  is  a  continuous  parameter  random  process 
(Ref.  22,  Sec.  3.7);  sample  values  of  8(t)  taken  at  some  time  intervals,  l/W,  are  describable  in 
terms  of  probability  by  a  joint  probability  density  function  p[8(l/W),  8(2/W), ....  8(n/W), . . .  ). 
Furthermore,  since  we  have  acknowledged  several  sources  of  distortion,  the  Central  Limit 
Theorem  suggests  that  the  real  and  imaginary  parts  of  8(t)  will  tend  toward  a  complex 
Gaussian  process.^ 


*v^  taken  m  positive  directed  oway  from  tfie  radar. 
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We  will  also  assume  that  the  process  is  stationary  in  time  for  a  Gaussian  process.  This 
requires  only  that  the  autocorrelation  of  samples  A(n/W)  and  d(m/W),  defined  as 

(where  the  bar  denotes  the  expectation  or  ensemble  average  of  the  process),  depend  only  upon 
the  difference  (n  -  m)/W,  and  not  upon  the  absolute  values  of  n  or  m. 

The  choice  of  the  distortion  model  of  the  form  of  Eq.  (8)  implies  that  the  distortion  vanishes 
with  the  signal.^  Furthermore,  bearing  in  mind  that  d(t)  is  complex  the  multiplicative  form  for 
the  distortion  proposed  indicates  that  the  distortion  amplitude  multiplies  the  undistorted  signal 
amplitude,  while  the  phases  add. 

The  signal  form  of  Eq.  (8)  will  therefore  be  used  throughout  this  report,  with  8(t)  taken  as 
a  stationary  complex  Gaussian  process. 

C.  The  Alternative  Hypotheses 

We  choose  to  frame  the  resolution  question  by  deciding  whether  or  not  a  second  signal  is 
present  in  the  side  lobe  region  of  a  first.  Formally,  we  are  going  to  choose  between  two  t^poth- 

jCt) 

eses  regarding  a  sample  of  received  complex  waveform,  CCt)  e  °  .  Specifically,  we  wish  to 
determine  whether  ((t)  is  from 

{(t)  =  u(t) +  <i(t,y^)  (10) 


or 


{(t)  =  U(t)  +  +  tit,-?  2} 


(11) 


where  v(t)  is  an  additive  noise  background.  The  noise  will  be  assumed  to  be  a  complex  Gaussian 
process  of  mean  zero  and  uniform  spectral  density  ("white"  noise)  of  2N^  watts/cycle,  for  posi¬ 
tive  frequencies  (if  one  physically  measured  the  noise  power  density  in  the  real  process  in  a 
band  the  result  would  be  N^^). 

Equation  (8)  appears  to  Ignore  the  fact  that  8^(t)  will  usually  be  related  to  8^;  the  larger  the 

target  return,  the  larger  (in  terms  of  probability  at  least)  the  distortion.  Consequently,  we  will 

relate  the  variance  of  8^  to  the  value  8^  (or,  if  8^  is  random,  to  its  variance).  Noting  that  the 

normalization  constraint  of  Eq.  (2)  relates  the  "no-error"  cross-section  8,  to  the  energy  E.  in 

th  ^  * 

the  real  part  of  the  return  from  the  i  target  by  the  factor 


=  E, 


(12) 


the  variance  in  8j(t)  will  be  taken  to  be  related  to  E^  by 

|«j(t)|^  =4Ej<r/  ,  (13) 

so  that  for  unit  cross  section,  the  variance  of  the  real  part  of  the  distortion  is  .  We  assume, 
for  the  present,  that  E^  is  a  known  constant  (in  later  sections.  E^  will  be  treated  as  a  random 
variable  and  we  will  replace  E^  by  E  in  the  equation).  The  mean  of  6^{t)  will  be  taken  as  zero. 

tMulHpolh  retvnw  con  be  cemMsrsd  a*  tsporale  tignoli  if  Ihe  polh  difference*  ore  large  enough  to  render  the 
above  CMumptien  lertowily  erroneou*. 
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While  we  will  investigate  the  effect  of  letting  be  a  random  variable,  it  will  be  assumed 
throughout  the  report  that  and  known  exactly.  This  assumption  is  not  really  as  restric¬ 

tive  as  it  might  be  at  first  assumed,  at  least  from  a  comparative  standpoint.  It  is  well  known 
that  if  the  problem  is  strictly  one  of  detection,  the  optimum  detector  for  processes  involving  un¬ 
known  target  parameters  can  be  theoretically  realized  as  the  weighted  sum  of  the  outputs  of  a 
bank  of  individual  filters.  Each  filter  the  bank  is  matched  to  a  different  combination  of  target 
parameters.  Thus,  the  heart  of  any  processor  for  signals  with  unknown  parameters  is  the  proc¬ 
essor  for  a  single  set  of  parameters. 

However,  an  apparently  important  distinction  lies  in  the  fact  that,  for  a  single  target  proc¬ 
essor,  we  need  only  test  all  possible  combinations  of  parameters,  while  for  even  pair-wise  reso¬ 
lution,  we  must  examine  all  possible  pairs  of  combinations.  Thus,  pair-wise  resolutions  neces¬ 
sitate  combining  the  outputs  of  approximately  n^/Z  filters  for  large  n,  while  the  signal  target 
processor  necessitates  combining  only  n.  The  price  we  pay  in  Increased  signal  energy  required 
for  this  increase  in  number  of  filters  appears  to  be  small,^  however,  if  we  can  draw  an  analogy 
with  analyses^^'^^  of  optimal  detection  of  one  of  M  possible  sirgle  signals.  For  example, 
Siebert's^^  analysis  of  the  Increase  in  E/N^  required  for  detecting  one  of  M  possible  orthogonal 
signals  indicates  that  the  increase  in  E/N^  required  with  M  is 

A(j^)«|lnM  .  (14) 

Bearing  in  mind  that  the  value  of  E/N^  which  is  required  to  detect  a  target  whose  positional 
parameters  are  known  exactly,  is  of  the  order  of  25,  and  that,  for  example.  In  20,000  «10,  it  is 
apparent  that  we  pay  little  penalty  for  dealing  with  one  of  M,  or  indeed  one  out  of  targets  for 
reasonable  M.  Consequently,  we  will  restrict  further  attention  to  the  "known  position"  case. 

In  summary,  we  will  concentrate  our  attention  on  the  pair-wise  resolution  of  targets  whose 
positions  are  known  exactly.  Only  the  distortion  of  the  signals  will  be  assumed  unknown  through¬ 
out  the  entire  report,  although  the  complex  target  cross  section  will  be  assumed  unknown  in  cer¬ 
tain  sections. 


D.  The  Effects  of  Signal  Distortion  on  the  AmUgulty  Function 


Before  considering  the  analysis  of  the  central  problem,  we  will  examine  the  effect  of  the 
signal  distortion  on  the  structure  of  the  ambiguity  function. 

For  distorted  signals  of  the  form  of  Eq.  (8),  it  is  useful  and  seems  appropriate  to  define^ 
the  ambiguity  function  as  the  cross  correlation  between  a  distorted  signal  and  an  undistorted 
reference  (thus  attaching  the  physical  significance  that  the  ambiguity  function  so  defined  is  the 
response  of  a  receiver  to  a  distorted  signal  with  parameters  when  the  receiver  is  matched 
to  an  undistorted  signal  with  parameters  72^- 


^(y^.  72)  -  ^12 


=  ^ ri»4  +  «j(t)) Pj(t) dt 


(15) 


t  Although  the  price  we  poy  in  hqrdwsre  complexity  certainly  may  net  be  miwII. 

t  it  will  be  OHwnsd  throughout  the  report  Ifwt  the  interval  0  to  T  it  long  enough  to  contain  all  tignal  tiructure 
of  intafott.  It  It  long  enough  to  juttify  the  replacement  of  the  llmilt  0  ^  T  on  intagrolt  by  -  oe  and  *  «e  with 
negligible  enor. 


10 


where  the  factor  normalizes  X  in  the  absence  of  errors  for  known  0^.  From  the  properties 
of  4j(t),  it  follows  that  the  mean  of  the  ambiguity  function  is  its  "no-error"  form 


J. 


12 


p^(t)  ^2^(1)  dt 


(16) 


and  its  variance  is  (recalling  that  the  real  part  of  X  is  of  concern,  and  that  the  variance  of  the 
real  part  is  one-half  the  variance  of  the  complex  X) 


(Xi2  X^2>  <^12  ^12>  ■ 


1^1 


V  ^  f^i  +  «i(t2)lp/(t,)  |i2(t^)  Pi(t2)  P2*(t2)  dti  dt2 


‘‘l  ^‘2<S>  »‘1<*2>  •‘2* <*2)  ‘**1  ^^2 


(17) 


which,  by  virtue  of  zero  mean  assumed  for  4^(t),  gives 
T  — : - 


2a. 


2  f*r  ^^1  6|(t2)]  ^  ,j, 

'  J  J  iTT^  **1  ^*1*  '^2*^1*  '^1**2^  >^2  *^2*  '^^l  '^*2 


1^1 


In  agreement  with  Eq.  (13)  for  |fl.(t)i  ,  we  will  set 


(18) 


<‘l’  *1<‘2>  =  -  ‘2> 

for  E.  known  and  constant,  where  p  is  the  normalized  [p(0)  =  1]  autocorrelation  function  of  the 

*  y 

errors.  Thus,  we  have,  recalling  that  E^  =  \0^\  /z, 

T 

"  ""6  ff  ^1  ^*1*  -  *2>  Pz  <*2>  *^*2  • 

0 

Let  us  attempt  to  attach  some  physical  significance  to  Eq.  (20).  Note  that  we  can  physically 
generate  the  ambiguity  function  by  the  use  of  a  linear  filter  of  impulse  response  h(t)  = 

P*(-t,  72^  ®  °  ,  as  indicated  in  Fig.  3,  if  the  interval  0  to  T  is  sufficiently  long  [a  finite  delay 

may  be  necessary  to  realize  h(t)|.  It  follows  from  the  properties  of  linear  filters  acting  on  a 
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3-41-5523 


LINEM  FILTER 
Ml) 


^(tI  ■  +  ,*(i,5j)/(t-T,J')dt 


■  VOI 


Fig.  3.  Linear  network  for  generating  X^2* 


Fig,  4(a).  Comparative  shapet  of  p(t)  and  p(t^  —  for  narrow-band 
distortion  of  coded  pulses. 


|3-«l-$5t«(ttT| 


Fig.  4(b).  Comparative  shapes  of  p(t)  and  p(t|  -  t-)  for  wide-bond 
distortion  of  coded  pulses  of  lorge-duratlon-bandwldth  product. 


12 


complex  Gaussian  process  of  this  type  (the  real  and  imaginary  parts  of  6{t)  are  independent,  with 
equal  variances)  that  the  output  of  the  filter  shown  is  Gaussian,  with  a  mean  value  which  is  the 
"no-error"  ambiguity  function,  given  by  Eq.  (16),  and  a  real  part  variance  given  by  Eq.  (20). 

This  result  makes  possible  a  useful,  practical  interpretation  of  a  rather  involved  integral 
which  will  appear  frequently  in  the  ensuing  analyses.  Since  X .  ,  is  a  complex  Gaussian  process, 
with  mean  value  ^^2*  the  amplitude  |X^2l>  is  distributed  according  to  a  "Rice"  probability 
density  function.^ 


where  Iq(x)  denotes  the  zero-order  Bessel  function  of  the  first  kind  of  imaginary  argument.  It 
then  follows  from  the  properties  of  this  function  that  at  points  where  the  "no-error"  ambiguity 
function  is  zero  (or  negligible  compared  to  the  error),  the  average  amplitude  of  the  ambiguity 
function  in  the  presence  of  errors  is  just  ^/2  Oy  Thus,  inspection  of  the  structure  of  X  from  a 
radar  allows  a  physical  estimate  of  Oy^ 

While  the  foregoing  fact  gives  physical  identity  to  the  integral  relation  of  Eq.  (20),  it  does 
not  lend  any  insight  into  the  relationship  between  and  Oy  For  specified  signals  and  autocorre¬ 
lation  functions,  the  integral  can,  at  least  numerically,  be  evaluated  with  precision.  Furthermore, 
the  integral  can  be  approximately  evaluated  for  limiting  cases  of  extent  in  time  of  p(t^  com¬ 
pared  to  the  correlation  time  of  ^(t)  (the  width,  in  time,  of  the  principal  maximum  of  the  ambiguity 
function,  »l/W). 

First,  if  p(t^  - 12)  is  essentially  constant  over  the  nonzero  interval  of  p(t),  as  indicated  in 
Fig.  4(a)  (corresponding  to  a  verj  narrow-band  distortion),  then 


indicating  that  <7^  «  in  the  ambiguity  function  side-lobe  region.  However,  p(t^  —  t^)  is  essen¬ 
tially  constant  over  the  entire  interval  of  nonzero  pit).  Thus,  such  error  represents  only  a  scal¬ 
ing  of  p(t)  and,  consequently,  a  scaling  of  X^^'  effect,  it  represents  a  component  of  unknown 
phase  and  amplitude  in  the  return.  Since  we  will  later  consider  signals  of  unknown  phase  and 
amplitude,  we  will  not  give  further  consideration  to  slowly  varying  p's. 

At  the  other  extreme  of  the  distortion  picture,  let  us  assume  the  distortion  to  be  quite  wide¬ 
band;  specifically  of  some  distortion  bandwidth  If  we  bear  in  mind  that  p{i^  -  is  the 

t  Alto  called  a  "modified  Rayleigh"  dentity  function  by  tome. 

t  The  problem  of  eitimating  ev  it  ttill  by  no  meotw  trivial .  A  major  obttacle  it  a  decition  about  the  thope  we 
thould  attribute  to  X]2.  For  many  practical  cotet,  the  tide  lobe  ttructure  of  the  ambiguity  function  will  be  olmott 
all  "error,"  in  the  tente  that  the  ttructure  will  bear  little  retemblance  to  the  "paper  detign"  thope.  Thit  tuggettt 
that  we  thould  tet  X)2  ^  ^  throughout  the  tide-lobe  region  at  o  practical  compromite,  and  rale|^  the  entire 
tide-lobe  ttructure  (averaged  over  tome  hopefully  repretentotive  templing  of  targeti)  to  tome  numerical  e^.  At 
the  other  extreme,  we  could  tet  atMe  five  minutet  et  each  hour  to  make  repeated  meaturenwnlt  of  the  omiguity 
function,  ilt  mean  and  ilt  variance  (on  tome  ihoughlfuily  provided  repreienlative  torgelt),  and  ute  Ihete  etti- 
nwtet  for  the  remainder  of  the  hour.  The  latter  procedura  would  undoubtedly  retult  in  a  much  tmaller  value  be-’ 
ing  ottigned  to  at  a  price. 
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Fourier  transform  of  some  spectral  density,  we  can  see  that  the  duration  over  which  p(t^  -  t^) 
is  large  will  be  of  order  l/W^.  Thus,  for  large  TW  coded  pulses  and  Wp  of  the  order  of  the 
modulation  bandwidth,  the  situation  will  be  as  indicated  in  Fig.  4(b).  A  useful  model  of  p(t^  -  t^) 
that  will  aid  in  visualization  (and  will  be  used  extensively  in  later  sections)  is 

p(tj  -t^)  =  sine  Wjj(t^  -  t^) 
where,  for  brevity,  we  follow  Woodward's®  notation 


This  form  of  p  corresponds,  of  course,  to  a  distortion  spectral  density  of  uniform  amplitude 
over  a  band  of  width  Wj^  cycles  centered  on  the  carrier  frequency,  and  zero  everywhere  else. 

As  p(tj  - 1^)  becomes  very  narrow,  it  acts  much  lilte  a  "delta  function"  of  area  l/Wjj,  and 
Eq.  (20)  becomes 

which,  for  the  form  of  p(t^  —  t^)  assumed,  for  T  very  large,  yields 


(21) 


From  the  normalization  constraint,  and  the  assumption  of  constant  signal  amplitude,  the 
integral  of  Eq.  (21)  is  expressible  as  the  ratio  T^/T^,  where  T^  is  the  time  overlap  of  Signals  1 
and  2  and  T^  is  the  total  time  the  radar  transmitter  is  on  for  each  signal  transmission.  Thus, 
for  broad-band  distortion  with  distortion  bandwidth  of  the  order  of  W^, 


<r 


2 

X 


^6 


Vo  ”^0 


(22) 


We  might  physically  expect  that  W^T  would  approximate  the  duration-bandwidth  product  for 
signals  of  this  type. 

It  is  thus  seen  that  for  a  fixed  value  of  <r^,  the  larger  the  duration-bandwidth  product,  the 
smaller  the  effect  on  the  ambiguity  function  side-lobe  structure  (a  corresponding  result  for  er¬ 
rors  in  antennas  and  their  patterns  is  well  Icnown).^^ 

For  a  given  distortion  variance,  the  variance  in  the  ambiguity  function  is  seen  to  increase 
linearly  with  the  degree  to  which  the  signals  overlap  in  time  for  coded-pulse  signals.  Since  it 
is  conunon  practice  to  concentrate  attention  on  the  case  of  appreciable  overlap  (T| «  T^),  cor¬ 
responding  to  the  region  of  the  ambiguity  function  close  to  the  main  lobe  in  time,  we  will  denote 
2  2 

the  value  of  for  this  case  by  ,  for  which  we  have 


Vo 


(23a) 


*We  might  tutpect,  however,  that  may  itself  tend  to  increase  with  increasing  TW. 
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and 


I 


(7 


2 

X 


(23b) 


For  pulse-train  signals,  a  form  of  p  that  seems  physically  justifiable  is  one  which  is  wide 
compared  to  the  width  of  an  Individual  pulse,  but  narrow  compared  to  the  interpulse  spacing,  as 
indicated  in  Fig.  5.  This  form  is  tantamount  to  assuming  essentially  constant  error  over  a  pulse, 
and  independent  errors  from  pulse  to  pulse.  The  constancy  of  error  over  the  pulse  represents 
a  reasonable  model  for  the  situation  in  which  the  individual  pulses  are  made  as  short  as  the  ra¬ 
dar  hardware  bandwidths  will  accommodate  (if  the  radar  bandwidths  are  everywhere  wide  com¬ 
pared  to  a  reciprocal  pulsewidth,  each  pulse  may  have  wide-band  distortion  of  the  type  previously 
treated,  but  this  "combination"  case  will  not  be  explicitly  treated  here). 

Under  these  assumptions,  t^  and  t^  of  Eq.  (20)  must  differ  by  less  than  one  pulsewidth  or  the 
integral  is  essentially  zero.  When  t^  - 1^  is  thus  small,  p(t^  - 1^)  can  be  ignored  in  the  integrand, 
and  letting  *2  "  *1  “  '*'*'**® 


where  d  denotes  the  pulsewidth  of  a  typical  pulse.  The  offset  in  the  factors  of  the  integrand  due 
to  T  just  makes  the  integrand  appear  to  be  the  convolution  of  two  pulse  trains,  K^(t|)  and  P^^t^), 
which  are  related  to  the  original  Pj(tj)  and 

lp;(t)|  =  IHi(t)l^  i  =  l,2 

(recall  that  a  binary  amplitude  is  assumed  throughout  the  train),  and  the  individual  pulses  of  the 
new  train  are  reduced  in  width  by  t,  as  indicated  in  Fig.  6.  The  value  of  lp(t)l^  is  l/T^  by  the 


Fig.  3.  AMumed  scale  of  the  error  correlation  function  p,  the  pulsewidth  d  and 
Interpulse  spacing  for  pulse  to  pulse  distortion  of  pulse  trains. 


ir 


l-SI-lill 


M(t,)- 


mi 


H  (♦,-▼! 


V*  vv*.  V 

»  vt*  0 

Mi. 


Fig.  6.  Visualization  of  products  of  the  form  p*  (t^)  (t^  -  t)  of  Eq.24. 
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normalization  constraint,  Eq.  (2),  where  =  Nd,  the  total  "on  time"  of  the  train.  A  first  in¬ 
tegral  (on  t^)  of  Eq.  (24)  produces  a  result  of  the  form 


(25) 


since  T ^d  =  N,  the  number  of  pulses  in  the  train. 

A  comparison  of  Eqs.  (25)  and  (22)  for  the  ambiguity  function  variance  for  pulse-train  signals 
and  continuous  signals  indicates  that  for  equal  we  have 


(Tj^  (pulse  train) 
<rj^(coded) 


—  if 


N 


The  product  of  distortion  bandwidth-time  duration  for  the  pulse  train  can  be  written  as  approxi¬ 
mately  ND/d  where  C  is  the  average  pulse  separation.  For  equal  products,  then 

(pulse  train)  k  T 

1  O 

Where  nonzero,  it  is  likely  to  be  of  the  order  of  l/N  or  more;  thus,  for  essentially  unity 

overlap. 


orj^(pulse  train) 
o’j^(coded) 


1- 

Nd 


Thus,  the  ratio  may  be  in  favor  of  either  waveform,  depending  upon  the  above  ratio. 

The  results  of  Eqs.  (22)  and  (25)  show  that  when  reasonable  care  in  alignment  has  been  taken 
and  small  target  distortions  are  present,  the  percentage  effect  of  the  signal  distortion  on  the 
ambiguity  function  will  be  negligible  in  the  main -lobe  region  of  the  ambiguity  function  (where 

«  1),  although  the  errors  may  completely  dominate  the  side-lobe  structure.  Thus,  the  error 
in  the  received  Signal  2  [<^(t,  of  (**))  will  have  negligible  effect  on  the  ability  to  detect  Sig¬ 
nal  2  in  the  presence  of  Signal  1 . 

Consequently,  since  we  will  concentrate  on  detecting  Signal  2,  it  will  be  assumed  that  Sig¬ 
nal  2  is  undistorted  for  the  remainder  of  the  report,  but  that  the  interfering  Signal  1  is  distorted 

j(t)  =  u(t)  +  [0^  +  «j(t))|ij(t)  t^^p^it)  .  (26) 


m.  THE  OPTIMUM  RESOLVER 

In  this  chapter,  we  will  derive  the  optimum  circuit  upon  which  to  base  a  decision  about  our 
two  alternative  hypotheses:  is  one  target  present  or  are  two?  We  will  analyze  the  performance 
of  this  resolver  in  terms  of  the  probability  that  we  will  correctly  say  Signal  2  present  when  it  is 
(the  probability  of  detection  Pq),  and  of  the  probability  that  we  will  say  Signal  2  present  wherf  it 
is  not  (probability  of  false  alarm  Pp). 
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A.  Deftnttlon  of  "Optimum" 

Before  attempting  to  derive  the  form  of  the  "optimum"  resolver,  the  sense  in  which  the  re¬ 
volver  is  to  be  (^timum  must  be  defined.  This  problem  has  been  discussed  at  length  by 
others,^' and  such  discussion  will  not  be  repeated  in  detail  here.  It  will  suffice,  instead, 
to  point  out  that,  as  one  might  expect,  the  most  we  can  ask  for  from  a  situation  which  is  basically 
probabilistic  is  further  probabilistic  information.  Consequently,  we  choose  to  consider  the  op¬ 
timum  receiver  as  one  which  computes  for  us  the  probability  that  one  of  the  hypotheses  is  true, 
given  the  available  data. 

The  hypothesis  we  wish  to  test  is  that  a  second  signal  present,  given  some  received 

waveform  ((t)  known  to  contain  some  randomly  distorted  signal,  which  we  will  denote  ^|(t).  In 
terms  of  probability,  we  wish  a  circuit  which  will  compute  the  a  posteriori  probability 
P[(('2(t)k(t),i/,^(t)]. 

If  ^^(t)  and  0^(t)  are  uniquely  defined  signals  [^2^^)  ~  02^^*' ^2^'  ^2  com¬ 

pletely  and  uniquely  specified],  we  can  make  straightforward  use  of  Bayes*  Rule^^'  and  write 
the  a  posteriori  probability  in  terms  of  the  a  priori  probabilities  (temporarily  suppressing  the 
time  dependence  in  the  interests  of  more  compact  notation)  as 

where  p[0|]  denotes  the  a  priori  probability  of  getting  the  particular,  unique  ^^(t)  of  concern,  and 
P^^l’  ^2^  denotes  the  joint  a  priori  prdMbility  of  getting  both  the  particular  and  il>^  of  concern. 

Since  any  monotonic  function  of  the  a  posteriori  probability  is  equally  (^itimum,  we  note  that 
Eq.  (27)  is  a  monotonlc  function  of  the  "likelihood  ratio"  L[(|i2(l))  given  by 

wz] "  -- itFTi;Tp(»;r  • 


A  useful  monotonically  related  function  is  In  L,  particularly  when  in  the  form 


In  **  i*' 


1  pftlll'l) 


since  this  form  separates  the  a  priori  probabilities  that  any  targets  and  are  present  into  a 
separate  factor  which  can  be  considered  as  part  of  the  threshold  setting  process  on  the  first  term. 
We  can  then  confine  our  attention  to  the  first  term,  which  we  will  denote  L'[^2(i))' 


L'[«2l 


=  In 


p[el»i.»2l 

p(tl^il 


(28) 


Our  primary  interest  will  be  directed  toward  the  situation  in  which  ^j(t)  and  ^2^^)  uniquely 

defined,  but  contain  random  compoiients  (the  random  distortion  process,  A(t)  of  Eq.  (12);  random 
complex  cross  sections,  p.).  Equation  (28)  and  its  development  remain  valid  for  this  case  if  we 
slightly  reinterpret  the  symbolism.  We  are  interested  in  the  a  posteriori  probability  of  some 
class  of  signals  given  a  sample  of  signal-plus-noise  C(t)  known  to  contain  one  of  an  en¬ 

semble  of  possible  <^|(t)'8.  If  the  statistics  of  the  random  parameters  of  are  known,  Eq.  (26) 
for  {(t)  allows  (at  least  in  principal)  the  straightforward  calculatiwi  of  p[i;k|j  or  ptckj,  <^2)' 
where  we  now  interpret  these  a  priori  probabilities  as  the  probability  density  functions  of  {(t), 
given  the  statistics  of  the  random  parts  of  ^^{\)  and/or 
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B.  The  Optimum  Resolver  -  All  Target  Parameters  Known 

In  order  to  first  obtain  some  kind  of  a  solution  to  the  problem  of  the  optimum  resolver  in  a 
comparatively  simple  context,  let  us  consider  a  highly  favorable  case  of  a  priori  knowledge; 

=  (/3j +«^(t)](i(t.yj)  , 

where  0^  7^  and  p(t)  are  known  exactly,  and 


^^(t)  = 

is  completely  known.^  Thus,  only  the  distortion  of  Signal  2  is  taken  as  random.  Then  £(t),  as 
given  by  Eq,  (26)  will  be  a  complex  Gaussian  random  process^  with  a  sample  probability  density 
function  of  the  form 


N  N 


2' 


..t 


-jTjexp  -i  ^  2  f«m-fm>^mnf^n“^nJ 

m=l  n=l 


where  is  the  i*^  sample  value  of  (;(t),  is  the  average  value  of  that  sample,  and  [A  ]  is  the 
covariance  matrix,  a  typical  element  of  which  is 


^kf  ■  ^^■k  “  ^k 

The  symbol  |  A  |  denotes  the  determinant  of  the  matrix  [A  ] . 

Assuming  the  samples  are  taken  at  regular  intervals,  spaced  some  At  apart,  we  have 
that  for  both  signals  present  (the  results  for  Signal  1  only  being  obtained  by  setting  0^  =  0), 

fk*^l»‘ik+^2f‘2k 

and 


(s;  -  fj)  (t,  -  f,)  =  (V,  +  )  (30b) 

where  the  subscripts  k  and  I  denote  samples  taken  at  t  =  kAt  and  lAt,  respectively.  For  white 
noise  of  spectral  density  2N^  watts/cycle  for  positive  frequencies,  and  zero  for  negative  fre¬ 
quencies,  the  autocorrelation  of  the  noise  is^ 

2NoUo[(k-f)  At)  (31) 

where  u^(t)  denotes  the  unit  impulse  function. 

We  then  have,  from  Eq.  (30b). 

^kf  =  2N^''oK‘‘-  At)  +4a/Ejp[(k-  I)  At)H*^Pj, 

where  it  is  recalled  from  Eq.  (19)  that 


*lk*ll  “  ■♦»/EjPl(k-f)  At) 


t  Recall  rtwt  it  wo*  shown  in  Sec. Il-B  that  the  effect  of  the  distortion  of  Signal  2  on  the  detection  of  Signal  2 
was  negligible  to  a  first  order. 
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where  p[(k  -  I)  At]  is  the  normalized  covariance  of  the  0(t)  process.  Writing  =  2E^/N^, 
we  have 


^ki  =  2^0  KK*'- -  O  At]  • 

The  problem  now  becomes  one  of  finding  the  inverse  of  the  A  matrix  A~^  such  that 


(32) 


N 

^  '^Im'^mn 
m=l 


8in»(i  —  n) 
»(f  -  n) 


(33) 


where  the  use  of  the  sinx/x  function  is  a  convenient  method  of  specifying  the  value  of  the  sum  to 
be  unity  for  t  =  n,  and  zero  otherwise. 

To  expedite  this  task,^  it  is  useful  to  make  the  notational  changes 

A^m  =  «{<^t.mAt] 


^mn  '  *‘[*"^*< 

and  note  that  (33)  can  now  be  written  as 

N 

I  «f£At.  mAt]  rfmAt.  nAt]  At  =  ^  "vflAt  -  ^tl 
m=l 


If  we  return  to  the  likelihood  ratio  of  Eq.  (28),  the  ratio  of  the  conditional  proltabilities  of  £  for 
Gaussian  processes  is 


P(fj.  •  •  • .  (-niff'll 


exp 


exp 


N  N 
S  E 
m,  n=i 
N 

S  S 


m,  n=l 


t^m-^m2»^mn2lgn-4l 
f^m  “  ^ml  J  ^mnl  ^^n  “  ^nl  • 


where  the  subscripts  1  and  2  on  the  statistics  of  £  refer  to  (1)  Signal  1  only,  and  (2)  Signals  1 
and  2  present. 

If  we  let  At  ->  0,  N  ->  <0,  such  that  NAt  =  T,  a  constant,  then  we  can  replace  the  double  sums 
by  integrals,  letting  mAt  =  t^,  nAt  =  t^,  to  yield 

p[e(t)|»,(t),»y(t)) 

p[t(t)U,(t)] 


exp 


exp 


"i  It *(t2)-t2<V  1*^*1 

0 _ 


T 


(35) 


where  r(t^,  t^)  satisfies  Eq.  (34)  in  the  limit.  To  find  this  limit,  we  note  that  the  area  under  the 
right  side  of  Eq.  (34)  is  unity,  regardless  of  At,  while  its  amplitude  approaches  infinity  at 


t  The  procedure  ihot  follows  is  based  upon  a  suggestion  by  W.M.  Slebert. 
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lAt  =  nAt,  and  zero  elsewhere  at  At  -*  0,  This  suggests  that  we  can  write  that,  in  the  limit, 
r(t2,  tj)  must  satisfy  the  integral  equation 

j*0  •  <36) 

For  A(t^,  t^),  we  need  only  make  notatlonal  changes  in  Eq.  (32)  to  establish  that 

«<‘l*  ‘e»  =  2N^f“o<*i  -  *2)  +  <n>  t‘l<‘2>I  • 

Let  us  assume  an  r(t2.  t^)  of  the  form 

'“<*2’  *3^  ”  2N~  f“o<*2  ~  V  ■*"  ^^<*2  “  *3^  ‘‘l  <*2^  til<*3H 
o 

where  A  is  assumed  to  be  a  constant  that  will  be  determined.  We  will  then  investigate  the  cir¬ 
cumstances  under  which  this  assumption  is  at  least  approximately  satisfied. 

Substitution  into  Eq.  (36)  shows  that  in  order  for  r(t2,  t^)  to  be  so  represented,  we  must  have 

<*1^  ►‘l<*3*  +  Ap(tj  -tj)  Pj  (t^)  p^ltj) 

+  ARj<r/p*(tj)pj(t3)  Ip,(t2)|^p(t^-t2)p(t2-t3)dt2  =  0  .  (37) 

The  value  of  the  Integral  will  be  considered  separately  for  two  different  classes  of  large 
TW  signals;  (1)  coded  pulses,  and  (2)  pulse  trains. 

For  coded  pulses,  assuming  the  distortion  to  be  wide-band  compared  to  the  reciprocal  of 
the  signal  duration,  the  integrand  is  of  the  nature  Indicated  in  Fig.  7.  It  is  api>arent  that  for  very 


Fig.  7,  The  itructwre  of  |p|(t^|^  p(t)  —  p(»2  —  13)  foe  wide-bond  dittortion 

of  o  pwlte  with  continuous  recrangulw  onvek^. 

large  duration-bandwidth  signals,  the  term  iP{(t2)|^  will  affect  the  value  of  the  integral  appre¬ 
ciably  only  when  t^  and/or  13  are  near  the  edges  of  the  pulse.  If  we  consider  the  pulse  to  last 
from  t  =  0  to  t  =  T^,  and  note  that  from  the  normalization  constraint  lHij(t)|^  = 
write  the  integral  as 

'*‘'2  -  *3>  «»»2  *  4;  Ph  -  4)  P<»2  -  h)  ‘>»2  ’ 


20 


p 


lll 

If  t2  and/or  is  outside  the  Interval  0  to  T^,  (t^)  and/or  111(13)  will  be  zero  and  the  equation 

is  satisfied,  regardless  of  the  integral.  Consequently,  if  the  distortion  is  truly  wide-band,  com¬ 
pared  to  l/T^,  the  limits  on  the  integral  can  go  from  -*»  to  •  with  little  error  (edge  effect),  and 
we  then  are  seeking  a  solution  to 

+  A]  p(t^  -  13)  +  T  ^  ^<*1  ■"  ~  ^^Z  =  °  ■ 

0  <  t^  <  T^,  0  <  tj  <  Tq  •  '  (38) 

It  is  obvious  that  this  equation  can  be  satisfied  if  p  is  of  such  a  form  that,  when  convolved  with 
itself,  the  result  is  again  p,  to  within  a  scale  factor.  Such  a  function  is  the  "sine"  function  in¬ 
troduced  in  Sec.  II-D; 

p(t^  -  tj)  =  sine  Wjj(t^  -  tj) 

This  choice  of  p(t^  -  t^)  implies  that  the  spectral  density  of  A(t)  is  uniform  over  some  distortion 
bandwidth  Wq,  a  not  unreasonable  model.  With  this  form  of  p,  Eq.  (38)  becomes 

u  o 

or,  from  Eq.  (23),  relating  and  <r^. 


A  = 


1  +  R.o  2 


Thus,  for  a  signal  with  a  rectangular  envelope  of  duration  T^  and  a  random  distortion  that 


has  a  flat  spectral  density  over  a  bandwidth  Wp  »  i/T^. 


r(tj,t2)»  jjj- 


R*a 


1  +  R,a^ 


(39) 


Combining  this  result  for  both  and  >'2(t|,t2)  of  Eq.  (33)  with  Eq.  (30a)  gives  the  form  of 

the  optimum  resolver  as  one  which  computes 


L«(*2<*)l  =  Re82 


I. 

» 


t(t)  P2  (t)  dt 


~  1  +  R  ff  ^  jy  ^**1*  **1  ^*1*  “  *2*  '*1*^2*  ‘‘2  ^^2*  ***2 


O  0 


(40) 


where  L"[^2(^)l  differs  fretn  L'[i^2(^))  ^9*  (^8)  only  by  a  factor  which  does  not  involve  the  re¬ 

ceived  data  of  C(t),  and  therefore  is  absorbed  into  the  question  of  threshold  settings.  Discussion 
of  the  significance  of  Eq.  (40)  will  be  postponed  until  the  optimum  resolver  for  pulse-train  signals 
is  derived. 
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Let  us  now  return  to  Eq.  (37)  and  pursue  a  solution  for  pulse-train  waveforms  of  the  type 
shown  in  Figs.  1(c)  and  1(d).  We  will  assume  that  the  Individual  pulses  are  not  coded,  and  that 
the  errors  are  highly  correlated  over  a  single  pulse,  but  essentially  uncorrelated  pulse  to  pulse. 
A  convenient  model  of  p(t^  - 12)  that  meets  these  requirements  is  a  Gaussian  shape 

p(t^-t2)  =expl--^-2-- 


as  indicated  in  Fig.  5,  where  a  is  assumed  greater  than  the  pulsewidth  d  but  much  smaller  than 
the  smallest  interpulse  period  D^.  Then  Eq.  (37)  is  approximately 


for  t^  and  t^  such  that  (t^)  ^  0. 

Under  the  above  assumptions,  the  Gaussians  involving  t^  and  t^  are  nearly  unity  when  t^ 
and  tj  differ  by  less  than  a  pulsewidth,  and  essentially  zero  for  t^  and  t^  separated  an  inter- 
pulse  period  or  greater.  Therefore,  for  t^  «  t^. 


+  A  +  ARj  ^ 


since  the  last  integral  is  approximately  the  ratio  of  the  area  under  a  typical  pulse  to  the  area 
under  the  entire  train  1/N.  Thus,  for  pulse  trains  of  equal-width  pulses 


A  = 


1  + 


N 


(41) 


Recalling  from  Sec.  II-D  that  we  can  establish  a  correspondence  between  the  number  of  pulses 
in  a  train  N  and  the  distortion-bandwidth  and  time-duration  product  comparison  of  this 

result  with  the  previous  value  of  A  for  coded  pulses  indicates  that  within  the  accuracy  of  the 
respective  approximations,  the  results  are  essentially  identical. 

With  this  identification,  Eq.  (40)  becomes  the  form  oT  the  optimum  detector  for  both  pulse- 
train  signals  and  those  with  continuous  rectangular  envelopes,  except  for  the  differing  forms  of 
p(t^  —  12)  assumed,  and  the  differing  constants  multiplying  the  double  integral. 

C.  A  Physical  Interpretation  of  the  Optimum  Resolver 

Inspection  of  Eq.  (40)  shows  that  the  resolver  consists  of  two  parts.  The  first  single  integral 
is  identifiable  as  the  output  at  time  zero  of  a  "matched"  filter  with  impulse  response 

*  a 

h(t)  =P2»‘2(-t)  e  ® 
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Identifying  the  double  integral  is  more  difficult.  However,  if  we  write  the  integral  aa 


^2  Jo  iC  ~  *^2  <^2>  ‘**2]  ‘*^1 


(where  we  revert  to  the  general  representation  for  the  normalized  autocorrelation  of  the  errors 
to  make  the  results  applicable  to  both  types  of  waveforms  under  consideration),  we  note  that  the 
term  in  brackets  is  just  some  function  of  t^  for  specified  signals  and  autocorrelation  function. 

For  example,  for  long  coded  pulses  and  wide-band  distortion,  the  "sine”  assumed  for 
p(t^  -  t^)  acts  approximately  as  a  delta  function  of  area  making  the  term  in  brackets 


^2  Jo  -t,)  p,(t,)  p*(t3)  dt, 

for  which  the  entire  double  integral  can  be  interpreted  as  the  output  at  time  zero  of  a  filter  with 
impulse  response 


h(t) 


D  O 


for  t  such  that  |pj(t)|^  ^  0 


0  ,  elsewhere 

Thus,  the  optimum  filter  for  this  case  is  one  with  an  impulse  response 


^2P2(-t) 


1  - 


1  + 


yrect  [u^(-t)l 


expljoj^t) 


where  we  define  the  function  "rect  Pilt)”  as  unity  during  the  duration  of  |p|(t)|,  and  as  zero 
elsewhere.  The  detector  is  thus  of  the  form  of  Fig.  8  for  this  type  of  signal  and  distortion. 

Fig.  8.  ^tlmum  rstoivsr  for  codsdpulis  lignoli 
with  wlds*oand  ditlortlon,  oil  peransisn  of  both 


l5-«l-55tll 


Cll|« 


•V 


LINEAfl  filter 


!(>|»TTE»wAToi?W  tlgnoliknownoxoetly.  Switch  In  Position  1  during 

nonzero.  Attenuation  by  factor 


t 


»/0+R,v^ ). 


During  the  portion  of  which  does  not  overlap,  we  convolve  with  ** 

would  if  Signal  1  were  never  expected.  During  the  overlap,  we  weight  the  filter  by  the  factor 


1  +  R,v  ‘ 


which  is  seen  to  have  a  value  of  unity  if  is  small.  Note  that  R^a  is  approximately  the 

ratio  of  the  mean-square  interference  of  Signal  1  at  Signal  2's  location  to  the  noise  power  density, 
if  the  overlap  is  very  large. 
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D.  The  Performance  of  the  Optimum  Resolver  -  All  Target  Parameters  Known 

In  this  section,  we  will  derive  performance  curves  for  the  optimum  resolver  of  Eq.  (40)  in 
terms  of  the  probability  of  detection  and  probability  of  false  alarm  Pp  for  a  Signal  2  in  the 
presence  of  Signal  1. 

We  note  from  Fig.  8  that  the  resolver  is  a  linear  device  (although  perhaps  time  varying)  so 
that  the  output  will  be  a  Gaussian  process,  as  is  the  input  process.  Thus,  we  need  compute  only 
the  mean  and  variance  of  the  output. 

We  assume  that  the  input  is  of  the  form  of  Eq.  (26).  Then,  from  Eq.  (40), 


L"(^(i2(t)l  =  Re 

+  0 


^2 

*  pT  * 

2  Jq  «i«t)p^(t)p2{t)  dt  + 


2 


i  +  r,.2 


T 

0 


T 

+  ^1^2  P(tl  -t2>  Pift2>  ^‘^*2> 

0 


T 

+  J  yj  «|(t|)|pj(tj)|^  p(tj  -  t^)  Pj(t2)  P2<t2>  <1*1  <1^2 
0 


T 

^  “l  '’‘*1  "  ^2>  “2  ‘‘^2 

0 


(42) 


when  p(tj  -  t^)  =  sine  Wjj(t^  - 1^).  If  Signal  2  is  present  the  mean  of  L"  is 


LJ  =  Re 


^1^2^12  *  1^2!' 


i  +  RjV  ‘ 


^1^2  “^2*  ►‘2*^2*  ''*1  **^2  ■'■  I'® 


2' 


(43) 


where  we  have  made  use  of  Eq.  (20)  in  the  last  term.  When  Signal  1  only  is  present,  the  mean 

A 

is  given  by  setting  P  ^  =  0  (but  note  that  P  ^  arises  from  the  form  of  the  resolver  and  is  not  set 
to  zero),  giving 
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"‘'l  ■  '2’  •‘*‘‘2’  •‘2‘'2’  "'1 


o  0 


(44) 


Thus,  the  presence  of  Signal  2  changes  the  mean  of  the  output  by 


AL"  =  \0 


,2'  ^<^0  Tq 

2  1  +  R^<t2 


(45a) 


where  we  have  let  p{t^  -  t^)  =  sine  Wjj(t^  —  t^)  and  assumed  "edge  effects"  negligible.  Reference 
to  Eqs.  (41)  and  (25)  indicates  that  the  form  of  (45a)  for  wide-band  distortion  of  pulse-train  sig¬ 
nals  is 


AL«  =  1^. 


‘  *  (^) 


(45b) 


1  +  R 


1  Ti 


12 


We  see  from  the  above  results  that,  for  either  Signal  2  present  or  absent,  in  notation  appro¬ 
priate  for  coded  pulses. 


L"  -  L»  =  Re  3 


tu(t) +aj(t)MLj(t)Jp2(t)  dt 

K  <r2 

■  7777^  ’ '  ** ‘'1^  <t2>  “^i 

1  n 


'O  0 


(46) 


from  which  the  variance  of  L  can  be  calculated,  noting  that  the  variance  of  the  real  part  of  the 
expression  is  half  the  variance  of  the  complex  process; 

T 

I  >i» 

<r/_  =  (L"-i7)2  = 


2 

-rA 


1  +  n^c' 


(f  ^ 


X  p(tj  -  t^)  p*(t2  -  tj)  ji^ttj)  '*^2  ‘**3 


*  (‘  ’  *■<)' 


r  T 


x»,(V  kJBj)  p(t,  -  tjl  -  t,l  bJ (t,l  Bjdjl  dt,  atj  dtj  t  r,b/  Jj’j’J’  |k,<1,)| 

0 

Xpjd^)  P2(t2)  p(tj  -t^)  P(t,  -tj)  P*(t3  -t^llpjitjll^pjd^)  dtj  dt^  dtj  dt^ 
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where  we  have  used  Eqs.  (19)  and  (20). 

For  the  long,  coded  pulse,  setting  p(tj  -  t^)  =  sine  Wpit^  -  t^),  assuming  wide-band  distor¬ 
tion,  and  neglecting  end  effects,  gives  after  some  algebra. 


=  1^: 


N 


1  + 


(47) 


Thus,  from  Eq.  (45)  and  Eq.  (47),  we  see  that  successful  detection  involves  choosing  correctly 
from  two  Gaussian  distributions  of  differing  means  and  equal  variances.  The  ratio  of  the  square 
of  the  difference  in  means  to  the  variance  is 


R 


•  *  (^) 


1  + 


(48) 


A  similar  analysis  for  pulse-train  signals  with  errors  correlated  during  the  pulse  and  decorre- 
lated  pulse  to  pulse,  results  in 

^1  -X, 


=R 
2  “2 


(49) 


12 


Note  that  both  results  can  be  written  in  the  form 


where  the  quantity  4  can  be  interpreted  as  an  "overlap  factor"  and  varies  both  with  the  degree 
of  time  overlap  T^/T  for  long  duration-coded  pulses,  and  with  the  "no-error"  ambiguity  func¬ 
tion  X^^  for  pulse  trains.  The  overlap  factor  is  plotted  in  Fig.  9  for  coded  waveforms. 

The  significance  of  this  parameter  can  be  appreciated  from  Pp  vs  Pp  considerations.  To 
compute  the  performance  curves,  we  need  only  note  that  for  T^  =  0,  we  have  the  classic  case 
for  detection  of  a  single  completely  known  signal,^’^*’^^’^^  and  the  Pp  -  Pp  curves  for  that  case 
are  applicable,*  with  only  the  values  relabeled.  Figure  10  is  a  set  of  such  curves  from  Manasse,^^ 
appropriately  relabeled. 

Note 'that,  in  Eqs,  (48)  and  (49),  we  can  set 


where  r  is  the  ratio  of  the  radar  cross  section  of  Target  1  to  Target  2.  Specializing  our  atten¬ 
tion  to  Eq.(48)  for  coded  waveforms,  note  that  if  -he  problem  is  essentially  the 


*Note  thgt,  for  a  single  known  signal,  the  inean_is  zero  when  no  signal  is  present,  whereas  Eq.(44)  indicates 
that  we  must  account  for  a  mean  s^ich  involves  Xin  in  the  threshold  setting  when  we  seek  a  secortd  target  over¬ 
lapped  by  a  first. 
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detection  of  R,  in  a  noise-only  background.  Consequently,  we  are  interested  here  in  the  case 
2  ^ 

»  1,  and  we  can  neglect  unity  in  the  denominator.  The  effect  of  the  interference  can  then 
be  assessed  by  setting  the  right  side  of  Eq.  (48)  equal  to  the  value  of  R^  t  desired  which  we  will 
denote  R^  (note  that  for  reliable  detection  R^  is  probably  on  the  order  of  50).  Solving  for  the 
value  R^  gives 


R 


2 


It  is  apparent  that  if  the  product  of  the  cross-section  ratio  and  the  "close-in"  ambiguity  function 
error  variance  is  much  greater  than  t/R^  (a  side  lobe  uncertainty  of  -30db  and  a  cross-section 
difference  of  20  db,  if  R^  =  50),  we  can  neglect  unity  in  the  numerator  and  write 


R 


2 


(50) 


E.  The  Optimum  Resolver  -  I%ase  and  Anq>litude  of  Signals  Unknown  - 
All  Other  Parameters  Known 

In  order  to  qualitatively  estimate  the  effect  of  our  assumption  of  complete  knowledge  of 
both  signals,  we  will  explore  in  this  section  the  effect  produced  by  lack  of  knowledge  of  either 
the  phase  or  the  amplitude  of  the  return  from  the  target;  first,  the  case  of  interfering  target 
cross  section  only  unknown,  and  then  both  cross  sections  unknown.  It  is  assumed  that  the  target 
cross  sections  are  samples  of  complex  Gaussian  populations.  The  two  cross  sections  of  interest, 
and  02<  further  assumed  to  be  Independent  and  of  mean  zero  and  to  allow  the  possibility 
that  they  arise  from  separate  populations,  we  will  assume  that  =  2E^q  and  -  ^E^q 

(note  that  we  will  hencefort^  distinguish  between  the  actual  energy  in  the  real  part  of  a  particular 
return  E^  and  its  statistical  average  E^^).  This  model  corresponds  physically  to  a  Rayleigh  am¬ 
plitude  distribution  and  a  uniform  phase  distribution. 

Since  the  similarity  of  results  for  both  coded-pulse  and  pulse-train  signals  has  been  demon¬ 
strated,  we  will  confine  attention  to  the  former  in  this  section. 

The  wide-band  distortion  process  will  be  assumed  to  be  a  complex  Gaussian  process  as  be¬ 
fore.  However,  in  order  to  solve  the  problem  with  any  degree  of  facility,  it  will  be  necessary 
to  assume  8^  and  8^(t)  to  be  independent,  except  that  we  may  choose  the  variance  of  the  8-process 
to  be  related  to  the  variance  of  in  a  manner  analogous  to  Eq.  (19): 


|«j(t)|^  =  4Ejj,(Tj^  .  (51) 

While  this  assumption  is  made  primarily  as  a  mathematical  convenience,  it  would  seem  to  repre¬ 
sent  a  "worst  case,"  in  that  any  known  correlation  could  be  capitalized  upon  to  our  advantage. 
Thus,  the  result  will  be  meaningful  -  at  least  as  a  limiting  case. 

The  simplest  approach  to  this  problem  is  to  consider  P  ^  fixed  at  first,  so  that  [compare  with 
Eqs.  (30a)  and  (30b) | 
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t(t)  = 
and 


«(tj  -  t2)  =  (t(t2)  -  {(tj)) 

=  j  u„(ti  -  t^)  +  ^  fl  +  Zcl  p(t^  -  t^)!  (tj)  ^i(t2)j  ,  (52) 

where  =  2Ej(j/N^.  If  we  assume  p(t^  - 1^)  of  the  form  sine  Wjj(tj  -  12)  and  ignore  end  effects, 
it  is  found  that 


*■<^1  ~  *2>  =  IN" 
o 


[l+RlO«T‘p(t,-t2)I  , . 

- ——2 - ,i4(t2) 

1  +  Rjo^, 


(53) 


If  we  now  form  the  logarithm  of  likelihood  ratio  for  a  fixed  P^,  we  find  that  the  terms  in¬ 
volving  the  received  data  are,  after  s<»ne  algebra, 

L"[*  I  =  Re  f  f  ^  pttt)  dt - --—2  t(t)  p!(t)  dt 

2 [Jo  1  +  Rj^a^  Jo 


-  ‘'1*  '’‘"l  ~  '‘1<‘2)  1^2  <*2> 


o  0 


(54) 


A  schematic  of  such  a  detector  is  indicated  in  Fig.  11.  Comparison  of  this  result  with  the 
optimum  detector  of  Eq.  (40)  (for  all  parameters  known)  indicates  that  for  targets  separated  so 


Fig.  1 1 .  Schematic  of  optimum  retolver  of  Eq.(54)  for  unknown  in  phese 
arid  amplitude. 
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that  results  approach  equality,  except  that  the  quantity  of  Gq.  (40)  is  replaced 

by  its  expected  value  R^^  in  Eq.  (54).  In  any  event,  the  resolver  of  Eq.  (54)  is  similar  in  struc¬ 
ture  to  that  of  Eq.  (40)  (Fig.  8).  except  for  one  additional  "box":  a  weighted  filter  matched  to  the 
undistorted  Signal  1. 

If  the  statistics  of  L"  are  examined,  straightforward  mathematics  will  prove  that 


1  -  1^12!^  +  Rio'^x^  (  °T  ) 


and 


‘'l= 


The  detection  parameter  of  interest  becomes 


(55) 


(56) 


(57) 


(^) 


'lO'^X 


(58) 


(compare  with  Eq.  (48) ) . 

Thus,  the  resolver  performance  is  given  by  the  curves  of  Fig.  9,  with  i  replaced  by 


V  = 


(59) 


Note  that  our  lack  of  knowledge  of  Signal  1  has  indeed  cost  us  something.  For  the  first  time, 
even  in  the  absence  of  errors,  we  must  pay  for  I X^  2 1  l>clng  large,  which  agrees  qualitatively 
with  Helatrom's  results^  (see  Sec.  I-A)  for  signals  with  unknown  parameters  but  known  X's.  Of 
course,  if  I X^  ^  I  ^  would  normally  be  true  if  the  targets  are  separated  the  width  ot  the 
main  lobe  of  the  ambiguity  function,  the  result  is  essentially  of  the  same  form  as  for  Signal  1 
known  completely,  except  that  R^  is  replaced  by  its  average  R^q. 

In  order  to  make  P  ^  appear  random,  we  need  only  multiply  the  antilog  of  Eq.  (54)  the 
a  priori  distribution  of  and  average  over  all 

r"  I  I 

®2ol  =  I  j  j  aE^o  j  ‘‘^a 
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The  result  is  the  same  as  the  corresponding  exercise  for  single-target  detection.®’ That 
is,  we  should  envel^e-detect  L"(4|i2)  of  Eq.  (54).  The  effect  on  detectability  is  well  known.®  For 
a  given  1^2 1  •  curves  of  Fig.  10  are  displaced^  to  form  those  of  Fig.  12  (see  Ref.  27)  where 
i'  is  given  by  Eq.  (59).  It  is  apparent  that  for  a  given  cross-section  magnitude,  the  additional 
price  which  we  must  pay  for  confessing  our  ignorance  of  exact  cross  section  in  the  resolver  de¬ 
sign  is  quite  small  at  levels  of  high  and  small  Pp. 

Thus,  it  is  apparent  that  for  large  overlap  (T,/T^  «  1 1,  resolution  with  high  probability  is 

*  °  2 

nearly  impossible  if  the  ambiguity  function  cannot  be  controlled  so  that  is  very  small.  How¬ 
ever,  for  this  optimum  resolver,  the  situation  improves  rapidly  to  the  point  that  for  as  much  as 
90  percent  overlap  ((T^  —  T^)/T^  =  0.1 ),  we  pay  only  a  10-db  penalty. 

A  comparison  of  Eqs.  (48)  and  (49)  indicates  that  if  we  can  realize  both  a  pulse  train  and  a 
long,  rectangular  envelope  signal  with  the  same  a^,  since  the  parameter  (1  ~^|2^  pulse 

train  can  be  made  to  drop  to  a  small  value  in  a  time  about  equal  to  the  width  of  one  pulse  of  the 
train,  the  pulse  train  might  provide  superior  resolution.  As  pointed  out  in  Sec.  II-D,  however, 
the  relative  value  of  <r^  that  can  be  expected  for  equal  care  in  hardware  design  and  alignment 
with  the  two  signal  types,  is  not  clear  at  this  time. 

IV.  THE  PERFORMANCE  OF  SOME  NCXfOPTIMUM  RESOLVERS 

This  study  would  be  incomplete  if  we  did  not  conclude  with  an  excursion  into  the  question  of 
"How  optimum  is  optimum?"  "Are  these  optimum  resolvers  sufficiently  superior  to  some  sim¬ 
pler,  nonoptimum  devices  to  be  worthy  of  practical  consideration?"  In  this  section,  we  will  at¬ 
tempt  to  at  least  partially  answer  the  question  with  some  examples. 

A.  The  Performance  of  an  Optimum  "No-Error”  Resolver  in  the  Presence  of  Error 

A  direction  which  suggests  itself  for  an  initial  venture  away  from  the  optimum  is  to  ignore 
the  error  in  the  signal  structure  in  the  design  of  the  receiver.  Specifically,  let  us  assume  that 
the  problem  is  that  of  the  previous  section:  resolution  of  targets  of  unknown  complex  cross 
sections,  but  all  other  parameters  known. 

Instead  of  constructing  exactly  the  filter  of  Eq.  (54)  (Fig.  11),  let  us  examine  the  effect  of 
assiuning  =  0,  and  of  constructing  a  filter  E^(^2l'  such  that 

Fl[^2(t)l=  |Fi(^2(t)ll  , 

where 

-^12t^i*<t)l  dt  (60) 

Then,  for 

t(t)  =  v(t)  +  (P,  +  «,(t)J  pj(t)  +  , 


t  Note  Ihot  a  different  Mt  of  curvet^  ^  it  generoled  if  we  inquire  into  Pq  vt  for  any  member  of  the  population 
of  cross  sections  characterized  by  the  assumed  distribution  for  ^2* 
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Fig.  12.  Pq  vs  Pp  for  aptimum  roiotvor  for  phoM  ottd  omplitudo  of  both  rolumt  unknown. 


Fig.  13.  Plot  of  ovorlop  factor  {  vt  ovorlop  in  rocoivod  lignoli  for  molvor  which 
noglacts  orror.  Tho  moon  intarfaranco-to-noito  ratio  I*  o  poromatar. 
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we  have 


r'T  * 

-Xj2j'J«i(t)|tii(t)|^dt+^2(l-  |Xi2l^)  - 


(61) 


where  it  is  noted  that  the  terms  Involving  iS^  add  to  zero.  The  statistics  of  F'  that  determine  the 
envelope  detector  output  are  the  mean  and  variance  for  a  fixed 


=  NohRl0<i7-l\2l'(*"«10"x  )1  • 

I  O  O  '  O' I 


Thus,  the  curves  of  Fig.  12  describe  the  vs  Pp  characteristics  of  this  "nonoptimum"  re¬ 
solver,  if  we  replace  the  4 '  of  Eq.  (59)  by  i  ",  given  by 


[1-  |x,2i^I^ 


*  *'l0‘^X  T  ~  1^2! 

o  o 


(62) 


V  +  «10^’ 


or,  for  large  overlap,  T^  «  T^,  and 


4"  « 


I-IX12I 


(63) 


The  overlap  parameter  for  this  resolver  is  plotted  in  Fig.  13  for  the  case  of 
Comparison  with  Fig.  9  for  the  overlap  parameter  for  the  c^timum  resolver  shows  that  for  large 
interference,  the  optimum  is  considerably  superior. 

We  can  quantitatively  compare  performance  by  asking  what  value  of  R2  is  required  to  achieve 
a  value  of  R2  4  "  =  Rg.  *nd  by  comparing  this  result  to  the  result  of  the  similar  example  in 
Sec.  III-D,  Eq.  (50).  Let  us  set 

Riq  =  rR^ 

and  assume  l^|2 1  ^  ^ 


«2<*  =  «o*«2 


2  ’^1 

*  +  “2*’%  T 
o  o 


Note  that  if  have 

o 
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1 


R  « 
O 


indicating  that,  no  matter  how  large  R^  may  be,  the  second  target  will  be  resolved  only  when 
the  overlap  is  small  enough  to  satisfy 


It  is  thus  apparent  that  with  this  form  of  resolver,  any  target  overlapped  appreciably  by  a 
target  for  which  the  product  ^  f  will  be  virtually  undetectable.  Thus,  we  can  say  that 

the  uncertainty  in  the  ambiguity  function  side-lobe  structure  sets  a  rather  definitive  limit  on  the 
difference  in  radar  cross  section  of  two  overlapping  targets  that  can  be  resolved  by  any  detector. 
Furthermore,  if  the  possibility  of  signal  distortion  is  ignored  in  constructing  the  resolver,  the 
limit  on  resolvability  is  quite  insensitive  to  the  radar  sensitivity  (signal-to-noise  ratio  achieved). 


B.  The  Resolution  Capability  of  the  Usual  "Matched”  Filter 

We  conclude  this  investigation  by  inquiring  about  the  degree  of  resolution  obtainable  if  one 
does  "nothing  special";  the  radar  processor  consists  only  of  a  filter  "matched*  to  which 

for  a  Gaussian  cross  section  is 

F2t*2<^>l= 

where 

£(t)p2(t)dt  .  (65) 

The  output  for  t(t)  =  u(t)  +  f/J^  +  «j(t))  ^^(t)  + 

T  T 

v(t)p2<t>  dt +^lX^2  +  «j(t)p,(t)p2(t)dt +/J2  .  (66) 

The  question  that  we  framed  was  what  happens  if  we  do  "nothing  special"  about  resolution. 
Obviously,  one  of  the  first  things  that  may  happen  is  that  the  response  of  the  filter  matched  to 
^2(1)  and  to  Signal  1  will  raise  the  false-alarm  rate  unless  the  threshold  is  appropriately  adjusted. 
To  see  this,  note  that  from  Eq.  (66).  the  variance  of  is 


indicating  that  the  threshold,  for  a  constant  false-alarm  rate,  must  be  raised  by  (in  terms  of 
voltage)  the  ratio 


T* 


T 


+  ff. 


o  o 


)1 


(67) 


(where  r  is  the  proper  threshold  setting  for  nonoverlapping  targets)  to  prevent  the  interference 
from  target  one  altering  Pp. 
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Since,  for  0^  fixed. 


we  see  that 


(F*)' 


R, 


1  +  R 


10 


from  which  the  potential  performance  of  the  simple  matched  filter  can  be  determined  from  the 
curves  of  Fig.  12,  replacing  by  4'",  given  by 


1  +  R 


'12' 


10 


T 

+  Rjo'^x  T 
o  o 


(68) 


It  is  apparent  that  if  I  ^  ^  1^12  ^  '^X  l^’^o  error 

in  X  predominates),  the  simple  filter  resolves  as  well  as  the  filter  which  ignores  only  the 
error  8.  (t)  examined  in  Sec.  IV-D.  However,  it  is  apparent  from  Eq.  (67)  that  the  requirement 
for  the  filter  mentioned  in  Sec.  IV-D  to  be  little  affected  by  ^^12^  «  1  if  the  targets 

overlap  appreciably;  for  the  matched  filter  it  is  necessary  that  ^'^X 

much  less  than  unity.  For  example,  for  »  T^,  we  can  write 


f  " 

^  (nonoptimum) 
^(matched) 


(1  -  \X 


12' 


or,  for  Rjo<7^ 


»  i. 


i  ” 

(nonoptimum) 

t  ni  * 

'(matched) 


(69) 


which  for  ay  «  l/z  and  |X.  ,|^  «  1,  as  we  would  expect,  is  clearly  greater  than  unity.  In 

I  £. 

regions  where  the  "no-error"  side  lobes  of  the  ambiguity  function  are  large  compared  to  the 
effect  of  errors,  the  filter  of  the  previous  section  is  markedly  superior.  This  result  is  not  un¬ 
expected,  since  the  major  difference  is  that  the  previous  filter  optimally  provides  for  the  "no¬ 
error"  ambiguity  function,  but  ignores  the  errors,  while  the  simple  matched  filter  ignores  both. 


V.  SUMMARY  AND  CONCLUSIONS 
A.  Summary  of  Results 

W''  have  explored  limitations  upon  the  ability  of  a  radar  to  detect  a  target  in  the  presence 
of  an  overlapping  target.  We  have  investigated  the  effect  of  the  inevitable  small  signal  distor¬ 
tions,  as  well  as  the  effect  of  lack  of  a  priori  knowledge  about  the  parameters  of  the  two  targets. 
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We  have  assiuned  that  the  signals  returned  from  the  target  are  randomly  distorted  by  a 
stationary  complex  Gaussian  process  8|<t)  so  that  the  distorted  signal  is  of  the  form 

(('^(t)  =  Idi  +  «i(t)]pj(t)  . 

where  /3|ti|(t)  is  the  undistorted  signal.  The  0(t)  process  is  assumed  to  have  mean  zero  and 
variance 


|«j{t)|^  =  , 

where  |^/2  is  the  energy  in  the  real  part  of  the  return  in  the  absence  of  error  (or  the 

expected  energy  if  the  cross  section  is  random). 

The  effect  of  this  distortion  is  to  contribute  an  uncertainty  to  the  shape  of  the  radar  am¬ 
biguity  function. 

We  showed  that  for  wide-band  distortion  of  coded  pulses  (see  Sec.  II-A  for  examples),  as¬ 
suming  a  flat  spectral  density  for  the  distortion,  the  variance  in  the  ambiguity  function  in  the 
region  near  the  main  lobe  was  related  to  the  signal  distortion  variance  by 


T  W„ 
o  D 


(23a) 


where  T^  is  the  pulse  duration  and  Wjj,  as  previously  pointed  out,  is  the  distortion  bandwidth. 
The  variance  in  the  ambiguity  function  varies  with  the  ratio  of  the  time  elapsed  from  main  lobe 
T^  to  the  duration  of  the  pulse  T^ 


(23b) 


where  we  can  interpret  T^/T^  for  the  resolution  problem  as  the  fraction  of  the  pulse  duration 
during  which  the  signals  overlap. 

We  found  that  for  pulse-train  signals,  and  for  errors  correlated  within  individual  pulses  but 
uncorrelated  pulse  to  pulse 


(T 


2 

\ 


(25) 


where  is  the  value  of  the  "no-error"  ambiguity  function. 

We  derived  the  optimum  resolver  and  calculated  its  performance.  We  found  that  a  convenient 
method  for  specifying  the  performance  of  any  of  the  resolvers  investigated  was  a  factor  4,  which 
we  designated  an  "overlap  factor."  since  it  defines  the  fractional  increase  in  signal-to-noise  ratio 
=  ^E^/N^,  required  to  compensate  for  the  presence  of  the  interfering  target. 

We  found  that,  for  the  optimum  resolver  as  specified  by  Eq.  (40)  when  both  signals  are  ex¬ 
actly  Itnown  (only  the  presence  of  the  second  is  in  doubt). 


4  = 


1  + 


(48) 
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for  coded  pulse  (see  Fig.  9)  and 


(49) 


for  pulse  trains.  We  found  that  the  value  of  required  for  Pjj  =  0.9,  Pp  =  10 
mately 

R  -  so 

2  "  T  -  T. 


-9 


was  approxi- 


(50) 


when  the  product  of  the  ratio  of  the  target  cross  sections  and  is  of  the  order  of  unity  or 

*0 

greater.  Also  we  found  that  for  the  target  cross  sections  unknown  and  assumed  to  be  complex 
random  variables,  so  that  the  mean  signal-to-noise  ratio  of  Target  1  R^q,  the  optimum  resolver 
is  an  envelope  detector,  operating  on  a  circuit  specified  by  Eq.  (54),  for  which 


1-  I5i 


i'  = 


I2I  +  (  °T  ‘) 


1  +  R,o<Tx 


(59) 


for  coded  pulses,  with  the  modification  for  pulse  trains  suggested  by  comparison  of  Eqs.  (48) 
and  (49).  The  curves  of  Fig.  12  are  applicable  for  R^. 

We  also  investigated  the  performance  of  two  nonoptimum  resolvers.  The  first  is  the  filter 
that  Eq.  (54)  indicates  would  be  optimum  in  the  absence  of  error  (we  assumed  in  the  filter  design 
that  =  0).  We  then  derived  the  performance  in  the  presence  of  errors,  and  found  that 


[i-  1^121^1^ 


^  *  ”l0*^X  T 
o  o 


12' 


(62) 


We  indicated  from  this  result  that  for  a  large  overlapping  target  resolution  can  be  achieved 
only  by  reduced  overlap,  and  not  by  increased  radar  sensitivity  (to  a  first  order). 

Secondly,  we  analyzed  the  performance  of  the  usual  "matched"  filter  (one  that  is  optimum 
for  nonoverlapping  signals).  We  pointed  out  that  the  presence  of  a  first  signal  will  cause  false 
alarms  if  the  threshold  is  not  prc^erly  varied,  according  to  Eq.  (67),  relative  to  the  ambiguity 
function  of  the  first  target.  If  the  threshold  is  properly  varied,  then 


**10*^X 


T 
o  o 


1  + 


(68) 


B.  Conclusions 

The  results  of  the  investigation,  notably  the  forms  of  Eqs.  (48),  (49)  and  (59),  indicate  that 
the  fundamental  limitations  on  the  resolvability  of  targets  are:  (1)  the  inevitable  unknowable  de¬ 
parture  of  the  shape  of  the  received  signal  from  that  expected,  and  (2)  our  lack  of  a  priori  knowl¬ 
edge  about  the  targets. 
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The  first  limitation  appears  to  be  the  more  restrictive  since  the  ambiguity  function  uncer¬ 
tainty  appears  in  the  performance  equations  always  in  the  form  of  R  .7  ^  ,  which  for  large  inter- 
fering  cross  sections  is  likely  to  be  considerably  greater  than  unity.'  For  the  distortion  model 
assumed,  the  uncertainty  in  signal  parameters,  at  least  as  indicated  in  the  difference  between 
Eqs.  (48)  and  (59)  above,  manifests  itself  only  in  the  factor  1  —  |\.  ,  { ^,  so  that  outside  the  main 
lobe  of  most  ambiguity  functions,  the  effect  is  certainly  negligible,*  regardless  of  the  cross 
section  of  the  interfering  target. 

The  hardware  cost  of  realizing  a  bank  of  (^timum  resolvers,  which  would  be  required  when 
signal  parameters  are  unknown,  would  probably  be  exorbitant.  Looking  toward  simpler,  non¬ 
optimum  resolvers,  an  interesting  comparison  can  be  obtained  by  comparing  the  performance  of 
the  three  resolvers  [Eqs.  (59),  (62)  and  (68)]  for  random  target  cross  sections,  for  the  situation 
in  which  interference  predominates  thermal  noise  *].  the  "no-error"  ambiguity 

function  magnitude  is  negligible  compared  to  unity.  The  quantity  which  seems  to  shed  most  light 
upon  the  performance  is  R2l.  which,  letting  the  target  cross-section  ratio  be  represented  by 
r  =  R^^R^  is  then  approximately 


for  the  optimum  resolver; 


R2«"  = 


r<r. 


o  o 


for  the  optimum  resolver  for  a  distortion-free  signal  and 


for  the  simple  matched  filter,  properly  thresholded.  Note  that  for  large  target  overlaps,  the 
quantity  la  independent  of  the  radar  sensitivity  (to  a  first  order,  at  least)  and  the  targets  are 
truly  unresolvable.  However,  the  optimum  detector  resolution  capability  increases  rapidly  with 
decreasing  overlap,  while  the  other  resolvers  improve  only  slowly  until  the  overlap  is  nearly  zero. 

The  superiority  of  performance  of  the  optimum  resolver  is  thus  clear.  The  choice  between 
the  two  nonoptimum  filters  revolves  about  the  relative  amplitudes  of  the  known  and  unknown  com¬ 
ponents  of  the  ambiguity  function.  It  would  appear  that  with  reasonable  s(H>histicati<m  in  design. 


t  Indeed,  If  it  im't,  thennal  noi«e  dominate*  Hie  detection  problem,  and  the  retolutlon  prabiom  i*  virtually 
nonexittent. 

t  Helihom's^  result*,  which  include  the  effect*  of  unknown  signal  orriyal  time,  would  seem  to  indicate  on 
increasingly  stronger  dependence  on  |1^|2 1  os  the  number  of  unknown  signal  parameters  increase*.  However,  his 
results  also  indicate  that  when  of  unknown  parameters  Is  negligible  in  on  optimum  r^ver. 
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it  should  be  possible  to  render  enough  over  most  of  the  ambiguity  function  space 

so  that  errors  will  predominate.  If  this  is  the  case,  the  simple  "matched  filter,"  properly 
thresholded,  would  appear  to  be  the  preferred  "nonoptimum"  detector  for  reasons  of  simplicity. 
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